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3 Abstract. The Gamma- Ray Observation of Winter THunderclouds (GROWTH) 

4 experiment, consisting of two radiation-detection subsystems, has been op- 
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erating since 2006 on the premises of Kashiwazaki-Kariwa nuclear power plant 
located at the coastal area of Japan Sea. By 2010 February, GROWTH de- 
tected 7 long-duration 7-rays emissions associated with winter thunderstorms. 
Of them, two events, obtained on 2007 December 13 and 2008 December 25, 
are reported. On both occasions, all inorganic scintillators (Nal, Csl, and BGO) 
of the two subsystems detected significant 7-ray signals lasting for >1 minute. 
Neither of these two events were associated with any lightning. In both cases, 
the 7-ray energy spectra extend to 10 MeV, suggesting that the detected 7- 
rays are produced by relativistic electrons via bremsstrahlung. Assuming that 
the initial photon spectrum at the source is expressed by a power-law func- 
tion, the observed photons can be interpreted as being radiated from a source 
located at a distance of 290 — 560 m for the 2007 event and 110 — 690 m 
for the 2008 one, both at 90% confidence level. Employing these photon spec- 
tra, the number of relativistic electrons is estimated as 10^ — 10^^. The es- 
timation generally agrees with those calculated based on the relativistic run- 
away electron avalanche model. A GROWTH photon spectrum, summed over 
3 individual events including the present two events and another reported 
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22 previously, has similar features including a cut-off energy, to an averaged spec- 

23 trum of terrestrial gamma-ray flashes. 
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1. Introduction 

24 Nonthermal X-ray and 7 ray emission, typically lasting for a few seconds to ~ 10 

25 minutes, has been observed from thunderstorm activity, with detectors on board an air- 

26 plane [McCarthy and Parks, 1985] and a ballon [Each et ai, 1996, 2000], high-mountain 

27 detectors [Suszcynsky et ai, 1996; Brunetti et ai, 2000; Chubenko et al, 2000; Alexeenko 

28 et ai, 2002; Muraki et ai, 2004; Toni et ai, 2009; Tsuchiya et ai, 2009; Chilmgarian et al, 

29 2010], and ground-based ones [Torii et ai, 2002; Tsuchiya et ai, 2007]. Interestingly, they 

30 do not appear to clearly coincide with lightning processes such as stepped leaders or re- 

31 turn strokes. In contrast, much shorter energetic radiation bursts, lasting only for tens 

32 of milliseconds or less, are often associated with lightning discharges. Though not neces- 

33 sarily homogeneous, they include terrestrial gamma-ray flashes (TGFs) [Fishman et ai, 

34 1994; Smith et al, 2005; Grefenstette et al, 2009; Briggs et al., 2010; Connaughton et al., 

35 2010; Marisaldi et ai, 2010a, b], natural lightning [Moore et ai, 2001; Dwyer et ai, 2005; 

36 Howard et al, 2008; Yoshida et al, 2008; Chubenko et al, 2009], and rocket-trigerred 

37 ones [Dwyer et al, 2003, 2004a, b]. 

38 In this way, it has recently become clear that apparently two types of radiation bursts 

39 with distinct duration are associated with thunderstorm activity. Although it is uncertain 

40 whether or not these two types have a common source mechanism, recent observations 

41 as well as theoretical works generally suggest that these bursts, especially short-duration 

42 ones, are produced by processes involving acceleration and multiplication of a background 

43 population of electrons. 



DRAFT 



January 12, 2013, 8:35pm 



DRAFT 



X - 6 TSUCHIYA ET AL.: GAMMA RAYS FROM WINTER THUNDERCLOUDS 

44 Various numerical kinetic calculations [Gurevich et al, 1992; Roussel-Dupre et al, 1994; 

45 Bell et al, 1995; Lehtinen et al, 1996; Gurehich et al., 1997; Milikh and Valdia, 1999; Gure- 

46 vich et al., 2007; Roussel-Dupre et al., 2008] and Monte Carlo simulations [Lehtinen et al., 

47 1999; Dwyer, 2003; Babich et al, 2005, 2007] commonly indicate that most of prompt non- 
48 thermal photons from lightning discharges are radiated, via bremsstrahlung, by relativistic 

49 electrons, which in turn are produced through mechanism involving relativistic runaway 

50 electron avalanche (RREA): some seed electrons, produced by e.g. cosmic rays, can be 

51 accelerated into relativistic regime if they can gain energies from the high electric fields in 

52 thunderclouds fast enough to overcome their total energy losses, due mainly to ionization. 

53 Then, they collide with air molecules and ionize them. Some of the faster newborn sec- 

54 ondary electrons are also accelerated to higher energies, hence increasing in their number. 

55 Finally, they will emit a detectable flux of nonthermal photons via bremsstrahlung. 

56 Early observations of long-duration bursts, though limited in number, measured X-ray 

57 fluxes in a few keV to a few hundred keV range, or 7-ray fluxes in MeV regions, suggesting 

58 that these prolonged emissions are also due to relativistic electrons [McCarthy and Parks, 

59 1985; Eack et al, 1996, 2000; Brunetti et al., 2000; Chubenko et al., 2000]. Several recent 

60 observations [Tsuchiya et al, 2007; Torii et al., 2009; Tsuchiya et al., 2009; Ghilingarian 

61 et al., 2010] have reinforced the suggestion, by detecting photon spectra extending clearly 

62 to 10 MeV or higher, and have given evidence that those long-duration 7 rays are also 

63 produced via bremsstrahlung. These results naturally lead to a view that long-duration 

64 events are also caused by relativistic runaway electrons. However, compared with short- 

65 duration ones, the nature of long-duration bursts have remained less understood, due 

66 primary to the lack of a sufficiently large sample. For example, it is still unclear how the 
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67 electron acceleration process keeps operating for such long durations. In addition, the 

68 relation between short-duration bursts and long-duration ones is unknown. 

69 Aiming at detections of radiation bursts from thunderstorm activity, we have been op- 

70 erating the Gamma-Ray Observations of Winter THunderclouds (GROWTH) experiment 

71 since 2006 December 20. In this paper, we report on successful GROWTH detections 

72 of two long-duration 7-ray bursts extending to 10 MeV. Using the acquired 7-ray data, 

73 the source distance, its spatial extent, and the number of relativistic electrons involved 

74 therein are estimated. Then, a 7-ray spectrum which sums up 3 GROWTH detections is 

75 compared with cumulative TGF spectra obtained by two independent space observations. 

76 Based on these results, we quantitatively discuss the production mechanism of prolonged 

77 7-ray bursts from winter thunderclouds. 

2. The GROWTH Experiment 

78 The GROWTH experiment, comprising two independent subsystems, has been oper- 

79 ating successfully at a roof of a building of Kashiwazaki-Kariwa nuclear power plant in 

80 Niigata Prefecture, Japan. Figure 1 shows the location of the plant, facing the Japan Sea, 

81 and the GROWTH experimental site therein. The geographical longitude, latitude, and 

82 altitude of the experimental site are 138°36 E, 37°26 N, and 40 m above sea level, respec- 

83 tively. This coastal area is frequently struck by strong thunderstorms in winter seasons. 

84 Actually, before the GROWTH experiment started working, radiation monitors (filled 

85 circles in Fig. 1), which are arranged at around 300 m — 400 m intervals in the plant, 

86 occasionally observed > 3 MeV intense radiation enhancements in winter seasons, which 

87 are difficult to ascribe to so-called radon washouts because these rainfall-related episodes 

88 would mainly cause increases at < 3 MeV energies [e.g. Yoshioka, 1992; Yamazaki et ai, 
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89 2002]. Each radiation monitor consists of a 05.1 cm x 5.1 cm Nal (Tl) scintillation counter, 

90 and a spherical ion chamber with a volume of ~ 14 L that contains Ar gas. The former 

91 covers the 50 keV — 3 MeV energy range, while the later operates in > 50 keV. However, 

92 the radiation monitors have too poor a time resolution of 30 sec, together with the too 

93 limited energy bands, to understand the nature of those phenomena. The GROWTH 

94 experiment is expected to provide much improved knowledge on these sporadic events. 

95 The pictures and drawings of the two subsystems are given in Enoto et al. [2007] and 

96 Tsuchiya et al. [2007]. One of them (Detector-A) uses two cylindrical Nal (Tl) scin- 

97 tillators (density = 3.67 gcm~^), having a diameter and a height of both 7.62 cm. In 

98 order to actively shield them from natural low-energy (< 3 MeV) environmental radiation 

99 (e.g. from ^°K), the Nal scintillators are individually surrounded by well-shaped BGO 

100 (Bi4Ge30i2; density = 7.1gcm~'^) scintillators, with the thickness on the side and bottom 

101 being 1.27 cm and 2.54 cm, respectively. The BGO scintillators geometrically shield the 

102 central Nal up to a solid angle of 2.4-7r str, or 0.6 x 4-7r. Thus, the Nal scintillators have 

103 a higher sensitivity toward the sky direction. The two central Nal scintillators and the 

104 BGO shields are operated over an energy range of 40 keV — 10 MeV. Output signals from 

105 photomultiplier tubes, attached to the Nal and BGO, are fed individually to a 12 bit 8ch 

106 VME-analog-to-digital converter [ADC (CP 1113A)] with a time resolution of 10 /isec, 

107 and is recorded on event-by-event basis. 

108 As another feature of Detector-A, a 0.5 cm thick plastic scintillator with an area of 

109 30.5 cm X 15.2 cm = 464 cm^ is placed above the two Nal scintillators, and operated with 
no a threshold energy of > 1 MeV. It has a high detection efficiency for charged particles, 
111 while it is almost transparent to photons due to its thinness and higher threshold. Utiliz- 
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112 ing this feature, we can separate charged particles from photons, and efficiently exclude 

113 background cosmic-ray muons, which typically deposit > 1 MeV energies, from events in 

114 the two Nal scintillators. Specifically, an event in either of the two Nal scintillators is 

115 judged as a charged particle if it give a simultaneous hit (with 10 fi sec) in the plastic 

116 scintillator. Thus, utilizing signals of the BGO and plastic scintillators both in anticoin- 

117 cidence, the central Nal scintillators effectively detect photons, generally arriving from a 

118 sky direction. 

119 Aiming at an independent radiation measurement, another subsystem (Detector-B) was 

120 installed ~ 10 m apart from Detector- A. It consists of spherical Nal (Tl) and Csl (Tl) 

121 scintillators (density = 4.51 gcm~'^), both with a diameter of 7.62 cm. The former op- 

122 erates in 40 keV —10 MeV, while the latter covers a higher energy range of 300 keV — 

123 80 MeV. Unlike Detector-A, these scintillators have omni-directional sensitivity because 

124 they have no shields such as the BGO or plastic scintillators. Output signals of two pho- 

125 tomultiplier tubes, attached to the Nal and Csl crystals, are sampled by a self-triggering 

126 electronics system with a 12 bit ADC (AD 574). These events are accumulated into an 

127 ADC histogram, which is recorded every 6 sec. 

128 Energy calibrations of Detector-A and Detector-B were carried out, using natural envi- 

129 ronmental 7-ray lines of ^i^pb (0.352 MeV), ^i^Bi (0.609 MeV), (1.46 MeV), and ^o^Tl 

130 (2.61 MeV). Then, especially for the Csl of Detector-B, cosmic-ray muons, giving energy 

131 deposits with its peak of around 35 MeV, were also utilized. Basically, these calibrations 

132 are performed and checked by an offiine analysis. 

133 In addition to those radiation detectors, the GROWTH system utilizes three optical 

134 sensors and an electric-field mill as environmental monitors. Each optical sensor consists 
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135 of a hand-made analog circuit, and a silicon photodiode (HAMAMATSU S1226-8BK) 

136 which is sensitive over a wavelength range of 320 nm — 1000 nm (with its peak at 750 

137 nm). They measure environmental visible light in coarsely different directions; sea side, 

138 zenith direction, and anti-sea side. The output signals are fed to a 12 bit VME-ADC, and 

139 recorded every 0.1 sec. The electric field mill is a commercial product (BOLTEK EFM- 

140 100). Its analog output is fed to a 12 bit ADC (AD 7892), and recored as electric-field 

141 strength between ±100kVm~^, with a resolution of 50 Vm~^. 

3. Results 

3.1. Count histories of the inorganic scintillators 

142 Figure 2 shows count histories of the 4 inorganic scintillators of Detector- A and B, 

143 obtained over 15:00 — 17:00 UT on 2007 December 13 which corresponds to local midnight 

144 (0:00 — 2:00 JST on 2007 December 14). For reference, typical background rates per 20 

145 sec, corresponding to the panels (a), (b), (c), and (d) of Fig. 2, are 22000, 1700, 2100, and 

146 1500, respectively. Similarly, Figure 3 gives those over 8:30 — 10:30 UT on 2008 December 

147 25 (17:30 — 19:30 JST on the same day, or local evening). On both these days, a strong 

148 low pressure system (with ~ 990 hPa on the ground) developed over Japan, causing 

149 thunderstorms at the coastal area of Japan Sea. A gradual count increase, followed by a 

150 gradual count decrease, generally shows that they are due mainly to radioactive radon and 

151 its decay products in rain,with their half-lives being 20 — 30 min. These effects originating 

152 from radionuclides are closely investigated by Suszcynsky et al. [1996] and Yamazaki et al. 

153 [2002]. 

154 Superimposed on such gradual count increases, a sharp count enhancement is found in all 

155 the inorganic scintillators at around 16:00 UT in Fig. 2, and at around 9:30 UT in Fig. 3. 
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156 Hereafter, we call the former and the latter events 071213 and 081225, respectively. These 

157 enhancements, both lasting for 70 — 80 sec, are quite different from the radon effects and 

158 from short radiation bursts associated with lightning discharges. Among those inorganic 

159 scintillators, BGO of Detector-A gave statistically the most significant burst detection on 

160 both occasions; SOcr for 071213 and 19a for 081225. This is because it has a higher density 

161 and a larger effective atomic number, and hence a higher stopping power, especially for 

162 X/7 rays, than the other inorganic scintillators used in our system. 

163 Figure 4 shows Nal count histories of 071213 in 3 energy bands from Detector-A and 

164 B, while Figure 5 represents those of 081225. For comparison with Detector-B, the data 

165 of Detector-A [panels (a), (b) and (c) of Fig. 4 and Fig. 5] are presented without the 

166 BGO or plastic anticoincidence. With a criterion that both the Nal and Csl scintillators 

167 of Detector-B simultaneously record 10 or higher counts per 12 sec in the 3 — 10 MeV 

168 energy band, we define burst periods of 071213 and 081225 as 84 sec, 15:59:29 — 16:00:53 

169 UT, and 72 sec, 9:28:29 — 9:29:37 UT, respectively. For reference, this energy band of 

170 either scintillator typically records ~ 4 — 5 events per 12 sec in quiescent periods; so the 

171 above criterion (again, not either but both scintillators have 10 or higher counts) means 

172 approximately > 3.2 — 4. Oct above the background. 

173 In order to estimate background levels of individual energy bands of Detector-A and B, 

174 we excluded data over the burst period (as defined above) and the adjacent 12-sec periods. 

175 The remaining data in the two lower-energy bands were fitted by a quadratic function 

176 (via evaluation), while those in the highest-energy band with a constant. Table 1 

177 summarizes the net count increases, obtained by subtracting interpolated background 

178 (dashed curves of Fig. 4 and Fig. 5) from the total counts in the burst period. Thus, 
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179 the burst detection is statistically significant in each of the three energy bands on both 

180 occasions. Table 1 also gives the observed photon number fluxes above the detectors 

181 using power-law spectra obtained later (Sec. 3. 5) and the detector responses of Detector 

182 B derived from a Monte Carlo simulation based on GEANT4 [Agostinelli et al, 2003]. 

183 Here, the MC simulation was evaluated with radionuclide sources of ^°Co and ^^''Cs. 

3.2. Arrival directions 

184 Shown in Fig. 4 and Fig. 5 are the Nal count rates of Detector A at > 3 MeV energies 

185 of 071213 and 081225, respectively, with [panel (d) in both figures] and without [panel 

186 (c)] anticoincidence. In both figures, the anticoincidence, which utilizes BGO and plastic 

187 signals in logical "OR", is seen to reduce the Nal background level (solid curves) to 

188 ~ 0.05 times that without anticoincidence. In contrast, the Nal-detected burst signal rate 

189 decreases due to the anticoincidence only to 0.25±0.03 and 0.31±0.06 times the raw rates, 

190 for 071213 and 081225, respectively. Thus, the burst photons survive the ant i- coincidence 

191 with 5 — 6 times higher efficiency than the background events. Similarly, the ratio of the 

192 > 40 keV Nal [Fig. 2 (b) and Fig. 3 (b)] to the > 40 keV BGO count rates [Fig. 2 (a) and 

193 Fig. 3 (a)], which is normally ~ 0.08 due mostly to environmental radioactivity coming 

194 from omni-directions, increased to 0.18 ± 0.02 for 071213 and 0.14 ± 0.02 for 081225. 

195 The above properties revealed by applying the anticoincidence are thought to reflect 

196 arrival directions of the burst signals. If they came mainly from horizontal or ground 

197 directions, the anticoincidence on/off ratio and the Nal/BGO ratio would both fall below 

198 their normal values, because, e.g., 40 keV or 3 MeV 7 rays horizontally entering Detector- 

199 A would be almost fully or partially (at least 30%) absorbed/scattered by BGO via 

200 photoelectric absorption and Compton scattering. Accordingly, we conclude that the 
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201 burst signals arrived from sky directions, not from horizontal or ground directions. A 

202 more quantitative study of arrival directions, employing Monte Carlo simulations, will be 

203 reported elsewhere. 

3.3. Burst components 

204 Figure 6 shows count histories of the plastic scintillator (> 1 MeV) and the envi- 

205 ronmental sensors. In coincidence with the apparent signals detected by the inorganic 

206 scintillators, the 0.5 cm thick plastic scintillator gave count increases in individual burst 

207 periods by A^pi = 160 ± 30 (5.3(t) for 071213 and 72 ± 18 (4a) for 081225 [top panels of 

208 Fig. 6]. Presumably these plastic signals are composed of 7 rays and charged particles, 

209 most likely electrons, which are either accelerated primaries or secondary ones produced 

210 by high-energy photons via Compton scattering. Below, we estimate how 7 rays and 

211 electrons contribute to A'pi, and estimate the electron flux above 1 MeV. 

212 First, a Monte Carlo simulation using GEANT4 predicts that the plastic scintillator has 

213 a low detection efficiency, 0.5 — 1%, for > 1 MeV 7 rays, while that for > 1 MeV electrons 

214 reaches 75 — 90%. Next, using power-law spectra obtained later (Sec. 3. 5) and the effective 

215 area of the Nal scintillator of Detector B yield 1 — 10 MeV photon number fluxes above the 

216 GROWTH system as (17.9±1. 8) X 10-2 cm'^s"^ for 071213, and (9.2±1.8) x lO'^ cm'V^ 

217 for 081213. Then, multiplying these fluxes by the GEANT4-derived detection efficiency 

218 for 7 rays of the plastic scintillator and the area of the plastic scintillator, 464 cm^, 7-ray 

219 produced counts contributing to A'^pi are estimated as A^^ = 70 ± 7 for 071213, and 31 ± 6 

220 for 081225. Finally, subtracting from A^pi, the contribution of electrons is obtained as 

221 90 ±30 for 071213, and 41 ± 19 for 081225. Although these numbers, when taken at their 

222 face values, imply a significant electron contribution to Api, here we conservatively regard 
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223 them as upper limits. Then, a 95% confidence level upper limit on the electron fiux above 

224 1 MeV of 071213 and 081225 is computed as 0.5 x lO'^ cm'^s"^ and 0.3 x lO'^ cm-^g-^, 

225 respectively. These upper limits are more than an order of magnitude lower than the 

226 1 — 10 MeV 7-ray fiuxes. Therefore, the observed burst signals arriving at the GROWTH 

227 system are inferred to be dominated by photons, rather than electrons. 

3.4. Comparison with signals from environmental sensors 

228 The visible-light sensor (middle panels of Fig. 6) recorded an extremely intense signal, 

229 lasting < 1 sec, at 15:57:50 UT and 9:25:09 UT, for 071213 and 081225, respectively. In 

230 coincidence with the recorded optical fiashes, the electric field rapidly changed its polarity 

231 from positive to negative (bottom panels of Fig. 6). These indicate that a lightning 

232 discharge occurred. However, their occurrence is well separated from the 7-ray bursts 

233 themselves, namely, 100 sec and 180 sec prior to the 071213 and 082125 commencements, 

234 respectively. Thus, we conclude that neither of the present two 7-ray bursts coincided 

235 with lightning discharges. 

236 Prior to the present work, Tsuchiya et al. [2007, 2009] have reported similar lack of 

237 coincidence between prolonged 7-ray bursts and lightning discharges. In Tsuchiya et al. 
23B [2007], a long-duration burst, lasting 40 sec, was detected by the GROWHT system 70 sec 

239 prior to lightning, while in Tsuchiya et al. [2009], no lightning discharges were measured 

240 over 5 minutes before or after a prolonged (~ 90 sec) burst detected at a high-mountain 

241 detector. These previously reported events, observed during thunderstorms, have also 

242 been considered to be associated with thunderclouds. 

243 In the same manner, we associate the present bursts with the thunderclouds, rather 

244 than to lightning discharges. Actually, rainfall-thunder observation data^, provided by 
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245 a laser observation system operated by Tokyo Electric Power Company, showed that 

246 thunderclouds approached the Kashiwazaki-Kariwa nuclear plant from the sea side on 

247 both occasions, and passed over it during the 10 minutes. 

3.5. Energy spectra 

248 Figure 7 and Figure 8 show background-subtracted GROWTH spectra, obtained in the 

249 burst periods of 071213 and 081225, respectively. In either case, we accumulated the data 

250 over the burst period, and subtracted background spectra which were averaged over 10 

251 minutes before and after the burst, although thunderstorms were ongoing during these 

252 time periods. This is to remove < 3 MeV line-7 rays induced mainly by radon decays, 

253 which increase the background level by up to twice. To examine how the background 

254 selection affects the final spectra (Figure 7 and Figure 8), we subtracted an alternative 

255 background spectrum averaged over 5 minutes before and after the burst. However, the 

256 background-subtracted spectra did not change by more than ±10% at < 1 MeV, or ±5% 

257 at > 1 MeV. These are almost negligible compared with the statistical errors. 

258 In both events, the background-subtracted spectra of Detector-A and Detector-B exhibit 

259 very hard continuum spectra, which clearly extend to 10 MeV. As shown in another 

260 GROWTH event 070106 reported previously [Tsuchiya et ai, 2007], and in high-mountain 

261 observations [Torii et al, 2009; Tsuchiya et al, 2009; Chilingarian et al, 2010], similar 

262 prolonged 7-ray emissions, extending to 10 MeV or higher, were observed, and have been 

263 thought to be produced via bremsstrahlung. Thus, the present high-energy 7 rays must 

264 also be produced via bremsstrahlung by electrons accelerated beyond 10 MeV. Given these 

265 results, the present two events, together with the previous ones, may be understood as 

266 manifestations of a common type of high-energy activity in thunderstorms. 
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As easily seen in Fig. 7 and Fig. 8, the obtained spectra, in particular those of Detector- 
B, flatten in 0.8 — 3 MeV, even though they are not corrected for the detector responses. 
One of the causes of this flat is the Compton scattering: since the Compton scattering cross 
section in the atmosphere increases as photon energy decreases toward 0.1 MeV, photons 
at low energies would experience stronger Compton degradation than higher-energy ones. 

3.6. Model fits 

Supposing that the burst 7 rays were produced in a source located at a certain distance 
and propagated through atmosphere to reach our detectors, we may deduce the initial 
photon spectrum at the source, and estimate the source distance, from the background- 
subtracted spectra. Since the Detector-A spectra are complicated due to the passive 
and active shielding effects by the BGO well, below we analyze the Detector-B spectra. 
According to numerical calculations [Roussel-Dupre et ai, 1994; Roussel-Dupre and Gure- 
vich, 1996; Lehtinen et ai, 1999; Babich et ai, 2007], an energy distribution function of 
runaway electrons, generated under the RREA mechanism, is expressed by a power-law 
function, or more precisely, an exponentially cut-off power-law. Consequently, we assume 
an initial photon number spectrum as 

/(ep) = ae-^exp {-eje,) (MeV-^sr-i). (1) 

Here, a and /3 are a normalization factor and a photon index, respectively, while ep and 
ec describe the emitted photon energy and a cut-off energy in MeV, respectively. While 
this equation represents an exponentially cut-off power law, it can also express a pure 
power- law by requiring ec -^^ 00. 
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276 Below, let us estimate the source distance d from our Detector-B data, as well as a, 

277 /3, and ec- In order to simulate the photon propagation in the atmosphere, we utilize 
27S EGS4 [Nelson et al, 1985] embedded in CORSIKA 6.500 [Heck et al, 1998]. In the 

279 CORSIKA simulation, the atmosphere consists of N2, O2, and Ar with the mole ratios 

280 of 78.1%, 21.0%, and 0.9%, respectively. The density of the atmosphere, divided into 5 

281 layers, depends exponentially on the altitude h, with a form of A + B exp {—h/C), with 

282 A, B, C being model parameters. For example, at h < A km, the model is specified as 

283 A = -186.6 gcm"^ B = 1222.7 gcm'^, and C = 9.94 km [Heck and Pzerog, 2009]. 

284 In addition, EGS4 can adequately treat electromagnetic processes in the relevant energy 

285 range of a few tens of keV to a few tens of MeV. 

286 Mono-energetic photon simulations were carried out for 33 incident energies from 50 

287 keV to 100 MeV. The energy interval is set to 10 keV for 50 keV - 90 keV, 100 keV for 

288 100 keV - 1 MeV, 1 MeV for 1 MeV - 10 MeV, and 10 MeV for 10 MeV - 100 MeV. 

289 For one mono-energetic photon simulation, one million photons were vertically injected 

290 to the atmosphere from a fixed source distance. In reality, 20 source distances from 20 m 

291 to 2000 m were applied for one mono-energetic simulation. Then, we saved the energy, 

292 angle, and species of all of photons and particles that arrive at the observatory level (40 

293 m above sea level). 

294 Figure 9 indicates three representative sets of simulated photon spectra, propagating 

295 over d = 300 m (36gcm-2), 1000 m (120 gcm'^), and 2000 m (220 gcm'^), with the 

296 numbers in parentheses giving air mass calculated by the above exponential formula. 

297 Punch-through photons, which suffer no interactions with air molecules, appear as a 

298 strong peak at the highest end of each photon spectrum, while scattered ones form a 
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299 continuum toward lower energies. As the distance increases, the punch-through photons 

300 and the scattered continuum are both strongly attenuated, in particular toward lower 

301 energies, due primarily to Compton scattering. For instance, the survival probability for 

302 10 (1) MeV punch-through photons to propagate over d = 2000 m is only 0.02 (10~^) 

303 times that over d = 300 m. Note that as discussed later in Sec. 4.3, a long-duration 

304 burst probably changes in the burst period its viewing angle relative to a beam axis of 

305 electrons accelerated in thunderclouds. Thus, the calculated photon spectra here would 

306 vary according to the changes, and hence they will be treated in this work as ones averaged 

307 over different viewing angles. 

308 Convolving the simulated photon spectra with the detector responses, we can obtain 

309 a model-predicted spectrum to be observed by the Nal and Csl scintillators. Finally, 

310 we convolve these model predictions with the assumed source photon spectrum, eq.(l), 

311 and fit the predictions simultaneously to the background-subtracted Nal and Csl spectra 

312 (right panels of Fig. 7 and Fig. 8). Then, the model parameters, such as a, (3 and ec, can 

313 be determined so as to minimize the fit x^- 

314 Shown in Figure 10 are three representative model fits to the spectra of 071213 and 

315 081225, assuming a power-law model. The choice of d of 300, 1000, and 2000 m in 

316 Fig. 10, respectively gave values as 48.4, 77.5, and 116 for 071213, and 40.9, 46.7, and 

317 51.4 for 081225. By changing d and repeating the fitting, we obtained curves as shown 

318 in Figure 11, together with the minima as 46.5 at d = 400 m for 071213, and 40.9 at 

319 d = 300 m for 081225. From these figure, we can constrain source distance. Table 2 and 

320 Table 3 summarize the best-fit parameters for 071213 and 081225, respectively, together 

321 with the constrained source distance. Also, low-energy parts (< 300 keV) of the two 
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322 spectra are found to play an important role to determine d. If we use the background 

323 averaged over the 5 min intervals (instead of 10 min), the source distances become 350 m 

324 for 071213, and 300 m for 081225. Thus, the distance is not affected significantly by the 

325 systematic background uncertainty. 

326 The Nal and Csl spectra have been explained, in either event, by a common set of model 

327 parameters, although the fits are not necessarily good enough. The cut-off energy was 

328 constrained to be rather high with relatively large errors. Thus, our data do not provide 

329 evidence for spectral cut-off in either event. In agreement with this, the two spectral 

330 models, a power law and an exponentially cut-off power law, gave similar goodness of fits 

331 in both events. Importantly, the source distance have been constrained with a reasonable 

332 accuracy. 

333 As another attempt, we tentatively fixed at 7 MeV, which is the expected average 

334 kinetic energy of runaway electrons (not of bremsstrahlung photons), and repeated the 

335 model fitting. Then, the fit became worse in both events (4th column of Table 2 and 

336 Table 3). Therefore, the initial photon spectrum is again inferred to extend beyond ~ 7 

337 MeV. 

4. Discussion 

4.1. Source heights 

338 Assuming a power-law function at the source, the 7-ray spectra of 071213 and 081225 

339 suggest that the sources are located at 290 — 560 m (35 — 67 gcm~^) and 120 — 690 m 

340 (14 — 82 gcm~^) above our system, respectively, both at 90% confidence level. In fact, 

341 these constraints are in good agreement with the known heights of winter thunderclouds in 

342 this area. Winter thunderclouds and winter lighting observed at the coastal area of Japan 
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343 Sea exhibit many features that have hardly been found from those in summer seasons 

344 and/or in other areas (e.g. Rakov and Uman [2005] and references therein). These include 

345 rather low altitudes of the development of these thunderclouds. Actually, Goto and Narita 

346 [1992] conducted video observations of winter lightning at the same Niigata Prefecture 

347 as our experimental site, and reported that the visible bases of winter thunderclouds are 

348 typically located at 200 m — 800 m above sea level. Also, a recent numerical calculation 

349 done by Babich et al. [2010] shows that another GROWTH event [Tsuchiya et al, 2007] 

350 may be produced at a source height of 0.5 — 2 km, and hence generally agrees with the 

351 present results. 

352 These height estimations provide an additional clue to the possible electron contribu- 

353 tions to the detected plastic signals (Sec. 3.3). As argued so far, electrons are considered 

354 to be accelerated in these thunderclouds to at least 10 MeV, probably a few tens of MeV. 

355 Since such electrons have a range of < 100 m at near the sea level, they would hardly reach 

356 our system, even if a range straggling is taken into account. Therefore, it is reasonable 

357 that the electron flux incident on our system, if any, was much lower than that of photons. 

358 Unlike the present sea-level observations, some high-mountain experiments, conducted 

359 at Mt. Norikura (2770 m) in Japan [Tsuchiya et al, 2009] and Mt. Aragatz (3250 m) 

360 in Armenia [Chilingarian et al, 2010], have detected primary electrons in long duration 

361 events (numbers in parentheses indicate altitudes of observatories) . Tsuchiya et al. [2009] 

362 estimated the source height as 60 — 130 m (90% confidence level), while Chilingarian et al. 

363 [2010] evaluated it as 100 — 150 m. These low source heights, which are comparable to or 

364 shorter than the expected electron range, can naturally explain their electron detections. 

4.2. Extent and motion of the 7-ray beams 
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365 Measuring electric-field structure of winter thunderclouds, Kitagawa and Michimoto 

366 [1994] revealed that tripole electrical structures, which consist of positive, negative and 

367 positive layers from top to bottom, appear at mature stages of winter thunderclouds. 

368 Then, they observed the tripole structures to last for < 10 minutes in early or late winter, 

369 while less than several minutes in midwinter. Since the present two events were observed 

370 in midwinter, the measured burst periods of 84 sec of 071213 and 72 sec of 081225 are 

371 consistent with their observations, if the burst durations represent the lifetime of electric 

372 fields. 

373 Figure 12 shows dose variations on 2007 December 13, measured by the nearest and the 

374 second nearest radiation monitors of the power plant (5 and 6 in Fig. 1, black and red lines 

375 in Fig. 12, respectively). The two monitors gave moderate dose increases for ~ 1 minute 

376 or less around the GROWTH event. By examining the GROWTH data of 071213 burst 

377 (crosses in Fig. 12), as well as dose rates of the second nearest monitor (red line) obtained 

378 for 15:59:30 — 16:01:00 UT and that of the nearest one (black line) obtained over 16:00:00 

379 — 16:01:30 UT, peak times of their enhancement can be evaluated as 15:59:48 (±6 sec) 

380 UT, 15:59:58 (±15 sec) UT, and 16:00:27 (±15 sec) UT. Thus, referring to the GROWTH 

381 data, the second nearest monitor increased in its dose rates with a small delay of 10 ± 16 

382 sec (or almost simultaneously), while the nearest one with a larger delay by 39 ± 16 sec. 

383 The two monitors are located at a distance of 500 — 600 m from the GROWTH system. 

384 For reference, data of the other two radiation monitors (4 and 7 in Fig. 1) exhibited no 

385 apparent increases (green and blue lines in Fig. 12). As for 081225, data of those radiation 

386 monitors were unavailable due to some data-storage problem. 
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387 These simultaneous and delayed detections by the two radiation monitors have two 

388 important implications. One is that the 7-ray emission from thunderclouds is likely to 

389 have illuminated a rather limited area, spreading over ~ 600 m on the ground. This kind 

390 of effect was also suggested by five radiation monitors (1—5 in Fig. 1), on the occasion 

391 of the other GROWTH event [Tsuchiya et ai, 2007], and another experiment conducted 

392 on the same coastal area [Torii et ai, 2002]. The other is that the 7-ray emitting region 

393 moved, presumably together with the thunderclouds. 

394 From data of Japan Meteorological Agency, it is found that south-west wind was blowing 

395 during 10 minutes including the burst period of 071213. Thus, the south-west direction 

396 can naturally explain the delay of the nearest monitor, if the 7-ray emitting region moved 

397 together with the thunderclouds. Then, the wind velocity was on average 360mmin~^, 

398 with the maximum of 720mmin~^. Projecting, to south-west axis, the distance between 

399 the GROWTH system and the nearest monitor, ~ 500 m, and dividing the projected 

400 distance, ~ 350 m, by its delay, 39 ± 16 sec, we obtain an average moving velocity of the 

401 emitting region as 540 ± 220 mmin"^. Thus, the estimated moving velocity is generally 

402 consistent with the wind velocity. 

403 Given above discussions, we may assume that the winter thunderclouds moved from the 

404 Japan sea in south-west side to the inland in north-east side. Then, a short-lived tripole 

405 structure appeared in a thundercloud, and accelerated ambient fast electrons toward the 

406 bottom positive layer. The accelerated electrons emitted 7 rays toward the ground, which 

407 the GROWTH system and the two radiation monitors detected when the beam passed over 

408 them. The differences in the statistical significance of detections between the GROWTH 

409 system and the radiation monitors may be due to different positions and effective viewing 
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410 angles relative to the accelerated electron beam axis in the thundercloud, and to the 

411 differences in their sensitivity. The 081225 event is considered to have occurred under 

412 similar conditions, because west winds, with almost the same velocity as in the 071213 

413 case, were blowing at that time. However, it is presently unclear whether the 7-ray 

414 emission ceased when the tripole structure disappeared, or when the 7-ray beam moved 

415 away from the GRWOTH system as the thunderclouds moved. 

4.3. The number of relativistic electrons in thunderclouds 

Using the initial photon energy spectrum /(cp) of Eq. (1) as quantified in Table 2 
and Table 3, we can estimate the number of relativistic electrons radiating the observed 
1 — 10 MeV 7 rays via bremsstrahlung, as 

iVe ~ T7 / / / n, sin^^g. (2) 

H Ji Ji Jo ri{Ke, ep, 0) 

416 Here, r]{Ke,ep,6) is the probability per 1 gcm~^ with which an electron with a kinetic 

417 energy emits a bremsstrahlung photon with an energy ep at an angle 6 with respect 

418 to the electron-beam axis [Koch and Motz, 1959], and H denotes the vertical length of 

419 the acceleration region. Since this H is unavailable from the present observations like in 

420 Tsuchiya et al. [2009], we assume either H = 300 m or 1000 m, corresponding to 35gcm~^ 

421 and llOgcm"^, respectively. These assumptions are based on intracloud observations of 

422 X rays using a ballon-born detector, which showed that a high electric field region, to 

423 produce a significant flux in 3 — 120 keV energy range, has a vertical extent of ~ 500 m, 

424 at altitudes of 3.7 — 4.2 km [Each et ai, 1996]. 

425 We further assumed that the electric-filed strength in the acceleration region is 300 

426 kVm~^, which is slightly higher than the threshold (at 1 atm) to cause the runaway elec- 
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427 tron avalanches. When 1 MeV electrons are accelerated from the top of this acceleration 

428 region to the bottom, they will gain energies of 15 MeV for H = 300 m, and 19 MeV for 

429 H = 1000 m. Therefore, the assumed electric-field strength, together with the assumed 

430 vertical length, is sufficient to produce 10 MeV photons via bremsstrahlung. 

431 To calculate eg. (2), we further need to specify 6; this is suggested to be relatively small, 

432 from the obtained photon spectra. In practice, the angle of prolonged 7-ray event may 

433 vary according to the motion of thunderclouds. Thus, we adopt 15° or 30° as O^^x- As 

434 listed in Table 4, these assumptions, together with eq. (1), give A'e = 10^ — 10^^. Similar 

435 estimations for other long-duration 7-ray bursts have given N^. = 10^ — 10^^ [Tsuchiya 

436 et ai, 2007, 2009; Chilingarian et ai, 2010]. Thus, long-duration 7-ray bursts appear to 

437 be emitted by a similar number of relativistic electrons. 

4.4. Relation between the bursts and the RREA mechanism 

A possible source of energetic seed electrons to cause the RREA can be attributable to 
secondary cosmic rays [Gurevich et ai, 1992]. The cosmic-ray flux above 1 MeV, at the 
presently relevant altitudes of <1 km, is Iq ~ 200m~^s~^ [Grieder, 2001]. Considering the 
measured burst periods of 70 — 80 sec and the 30 — 40 sec delay of one radiation monitor 
from the 071213 event, we presume that the acceleration region has a horizontal length 
of L ~ 600 m at most, as judged from the extent of the 7-ray beam of 071213. Possibly, 
an actual extent of the acceleration region in thunderclouds would be shorter than this 
600 m, because the 7-ray beam would diverge due to multiple scatterings of the emitting 
electrons and Compton scatterings of the emitted 7 rays. We may also consider that an 
acceleration region is sustained in thunderclouds at least for 100 sec. Accordingly, the 
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number of cosmic rays Sq, entering the acceleration region, is described as 

^0 = 7.2 X 10^ X (L/600 m)^ x At/100 sec. (3) 

Based on the RREA mechanism, the total number of relativistic electrons at the end of 
an acceleration region, Wrrea, is estimated as 

iVRREA = ^oexp(5), 6= —. (4) 

JO A 

The length parameter A is given as 

7300 kV 
-C/ — (276 kV m 

where E is the electric-field strength in kV m~^ and n denotes the air density relative to 
that at 1 atm. This formula is valid for 300 — 3000 kVm^^ [Dwyer, 2003]. Assuming 
E = 300 kV m~^ gives A ~ 300 m at P = 1 atm. In practice, E may be somewhat lower 
than 300 kVm~^, because P during thunderstorms would be usually lower than 1 atm 
due to lower pressure system, and hence gives n < 1. Since a uniform field, S = H/X, 
gives A^'rrea = Soexp {H/ X), the factor r] = exp{H/X) is regarded as the avalanche 
multiplication factor, and becomes 3 and 30 for H = 300 m and 1000 m, respectively. As 
a result, we obtain 

A^RREA = 2.2 x 10^^ X (L/600 m)^ x At/100 sec x r//30. (6) 

We thus obtain A'rrea = 10^° — 10^\ which agrees generally with the derived = 
10^ — 10^^. Thus, the standard RREA process can explain at least the present two 
prolonged bursts. 

In the above estimation, we assumed an electric field is slightly higher than the RREA 
threshold. However, a weaker field below this threshold might suffice to produce prolonged 
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443 7-ray emission. In reality, a 30 — 120 keV x-ray flux continuously increased while an 

444 electric field is lower by 30% — 60% than the RREA threshold [Each et ai, 1996]. This 

445 quasi-static moderate-level field might be accomplished by e.g. a charging mechanism of 

446 thunderclouds. 

4.5. Comparisons with TGFs 

447 The derived (sec. 4. 4) is more than five orders of magnitude lower than the number of 

448 relativistic electrons expected from TGF observations, e.g. 10^^ — 10^^ [Dwyer and Smith, 

449 2005]. This huge number of relativistic electrons in TGFs may be generated by relativistic 

450 feedback mechanism, involving positrons and X rays propagating in the opposite direction 

451 to runaway electrons [Dwyer, 2007, 2008]. Since the estimated A'c of the present bursts 

452 is in generally agreement with A^rrea expected from the simple RREA mechanism, we 

453 conclude that at least the present two events do not require an intense feedback process. 

454 In order to better characterize 7-ray spectra of long-duration events, we stacked count 

455 spectra over three bursts, namely, the present two ones and 070106 [Tsuchiya et ai, 

456 2007]. Figure 13 compares the summed GROWTH spectrum with averaged TGF ones 

457 obtained by two independent satellites; one sums 289 events measured by Reuven Ramaty 

458 High Energy Solar Spectroscopic Imager (RHESSI) [Dwyer and Smith, 2005], while the 

459 other averages over 34 events observed by the Astrorivelatore Gamma a Immagini Leggero 

460 (AGILE) satellite [Marisaldi et ai, 2010a]. Thus, the composite GROWTH spectrum is 

461 rather similar in shape to those from TGFs, although three spectra all include detector 

462 responses. This resemblance is consistent with our basic standpoint [Tsuchiya et ai, 

463 2007, 2009] that the long-duration 7 rays are emitted by the same bremsstrahlung process 

464 as TGFs. 
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465 On a close comparison, the TGF spectra, especially the AGILE one, appear to have 

466 a higher cut-off energy than the GROWTH spectrum. This may be attributable to a 

467 difference in the electric potential operating in an acceleration region [Dwyer and Smith, 

468 2005]. Thus, electrons accelerated in a much lower atmospheric density at the production 

469 sites of TGFs, 15 — 40 km [Dwyer and Smith, 2005; Carlson et ai, 2007; 0stgaard et al, 

470 2008], would propagate through a longer distance, which gives a higher electric potential, 

471 and gain higher energies because of a smaller ionization loss per unit length. 

5. Summary 

472 The GROWTH experiment observed two long-duration 7-ray emissions from winter 

473 thunderstorms on 2007 December 13 and 2008 December 25. The photon spectra ob- 

474 tained in both events clearly extends to 10 MeV, and are consistent with a scenario 

475 that accelerated electrons produce, via bremsstrahlung, the observed 7 rays. Adopting a 

476 power-law function as the initial photon spectrum at the source, we have constrained the 

477 source distance as 290 - 560 m for 071213 and 110 - 690 m for 081225, both at 90% con- 

478 fidence level. These constraints agree with visible-light observations, which show that the 

479 bottom of winter thunderclouds is usually located at 200 — 800 m above sea level [Goto 

480 and Narita, 1992]. We have shown a possibility that the observed 7-ray beams move with 

481 winter thunderclouds, and spread over ~ 600 m. 

482 We estimated the number of relativistic electrons to cause the present prolonged 7-ray 

483 emissions as 10^ — 10^^. These are in general agreement with those expected from the stan- 

484 dard RREA mechanism triggered by secondary cosmic rays. The cumulative GROWTH 

485 spectrum, summed over the present two ones and another GROWTH event [Tsuchiya 
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486 et ai, 2007], was found to be similar in basic spectral features with the averaged TGF 

487 spectra [Dwyer and Smith, 2005; Marisaldi et ai, 2010a]. 
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Figure 1. The location of tlie Kasliiwazaki-Kariwa nuclear power plant (left) and its bird 
view (right). A filled square in the right panel represents the GROWTH experimental site, while 
9 filled circles show locations of radiation monitors. Each original image is taken from Google 
Map. 
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Figure 2. Count rates per 20 sec of the 4 inorganic scintillators over 15:00 — 17:00 UT on 
2007 December 13. Panels (a) and (b) show the > 40 keV count rates from the BGO and Nal 
scintillators without the anticoincidence of Detector-A, respectively, while panels (c) and (d) 
represent those of Nal (> 40 keV) and Csl (> 300 keV) scintillators of Detector-B, respectively. 
Horizontal axis shows universal time. Error bars are statistical la. The gaps in panels (a) and 
(b) are due to regular interruptions of data acquisition of Detector-A every hour. 
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Figure 3. The same as Fig. 2, but over 8:30 - 10:30 UT on 2008 December 25. 
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Figure 4. Count histories per 12 sec of 071213 in three energy bands from Detector-A (left; 
summed over the two Nal) and Detector-B (right; Nal), obtained for 15:55 — 16:05 UT. Panels 
(a), (b) and (c) correspond to 0.04 — 0.3 MeV, 0.3 — 3 MeV and > 3 MeV energy bands without 
anticoincidence, respectively. Panel (d) indicates the > 3 MeV energy band of the Nal with 
anticoincidence. Panels (e), (f) and (g) are the same as Panels (a), (b) and (c), respectively, 
but for Detector-B. Solid curves outside the burst period show the estimated background level 
(see text), while dashed ones denote the interpolated background level over the burst period. 
Abscissa represents universal time. Each error bar is statistical la. 
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Figure 5. The same as Fig. 4, but for 081225, obtained over 9:24 - 9:34 UT. 
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Figure 6. The count-rate histories of the plastic scintillator of Detector A and the environmen- 
tal sensors. Left panels show 071213, obtained over 15:55 — 16:05 UT, while right ones represent 
081225, obtained for 9:24 — 9:34 UT. Top, middle and bottom panels in both sides represent 
a > 1 MeV counts every 12 sec from the plastic scintillator, 1-sec optical data variations, and 
1-sec electric field ones, respectively. All abscissa are universal time. Vertical lines in top panels 
represent the burst periods. 
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Figure 7. Background-subtracted spectra of Detector-A (left) and B (right) of 071213. Black 
and gray points in the left panel indicate the Nal data without and with anticoincidence, respec- 
tively, while those in the right panel show the Nal and Csl scintillators, respectively. All error 
bars quoted are statistical la. Arrows, showing 95% confidence level upper limits, are draw- 
ing when statistical significance of a data point is lower than la. The horizontal and vertical 
axes show the photon energy in MeV and count per unit energy interval, respectively. Detector 
responses have not been removed. 
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Figure 8. The same as Fig. 7, but for the 081225 event. 
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Figure 9. Photon spectra at the observatory derived from Monte Carlo simulations. Three 
panels correspond to = 300 m (left), 1000 m (middle), and 2000 m (right). Different colors 
denote incident photon energies, Eq = 0.3 (blue), 1 (green), 3 (red), and 10 MeV (black) in 
all panels. Abscissa shows the photon energy at the observatory, while ordinate represents 
probability density function. 
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Figure 10. The photon spectra observed by the NaT (filled circles) and Csl (open circles) 
scintillators of Detector-B, compared with calculations for assumed source distances of 300 m 
(black), 1000 m (red), and 2000 m (green). Left panel shows 071213, while right one indicates 
081225. For clarity, the Csl data and the corresponding model spectra are multiplied by 0.1. The 
horizontal and vertical axes show the photon energy in MeV and counts in each bin, respectively. 
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Figure 11. The values of A^^ = ~ Xmin? plotted as a function of the assumed source 
distances (black circles). Left and right panels show 071213 and 081225, respectively. Black 
curves show smoothing lines. Horizontal dashed lines from bottom to top correspond to 68%, 
90%, and 99% confidence level. 
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Figure 12. Radiation dose rates per 30 sec (left ordinate) obtained by ion chambers of 



radiation monitors over 15:55 — 16:05 UT on 2007 December 13. Different colors specify the 



radiation monitors numbered as 4 (green), 5 (black), 6 (red), and 7 (blue) in Fig. 1. Superposed 



on them, crosses show the > 40 keV count history per 12 sec of the Nal detector of Detector-B 



(right ordinate). Vertical dashed lines represent the defined start and end time of the 071213 



event. 





-1 

■L 






p-, [—1+ ^^L= 





DRAFT 



January 12, 2013, 8:35pm 



DRAFT 



TSUCHIYA ET AL.: GAMMA RAYS FROM WINTER THUNDERCLOUDS 



X - 45 



>10 f 
« 10 



,2 ^ 



§10-^ 
o 

10-2 



I • This work ^ 



10"^ r A AGILE 
A RHESSI 



10"^ 1 10 

Photon energy (MeV) 

Figure 13. Cumulative spectra of tlie tliree GROWTH events from Detector-B [Nal and 
Csl (filled and open circles, respectively)], compared with summed TGF spectra by RHESSI 
and AGILE. The latter two spectra are adopted from Fig. 2 of Dwyer and Smith [2005] and 
Fig. 5 of Marisaldi et al. [2010a]. For clarity, the RHESSI and AGILE spectra are multiplied 
by 1 X 10~^ and 1 x 10~^, respectively. The vertical and horizontal axes represent counts in 
MeV~^s~^ and photon energy in MeV, respectively. Errors assigned to the GROWTH and AGILE 
data are statistical ones, while those for RHESSI include systematic uncertainty of background 
estimation [Dwyer and Smith, 2005]. 
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Table 1. The count enhancements and the corresponding photon number flux of 071213 and 
081225. 

071213 081225 

(MeV) Det. Det. flux=(xlO-^cm-^s'^) Det.A^ Det. B'' flux^(x IQ-^cm-^s"^) 

0.04 -0.3 530 ±70 1120 ± 110 20.2 ± 1.9 300 ± 60 390 ± 90 6.1 ± 1.7 

0.3 - 3 900 ±80 1660 ± 140 22.2 ± 2.3 400 ± 60 920 ±110 7.6 ±2.1 

3-10 410 ±30 370 ±20 10.5 ±1.0 180 ± 20 178 ± 16 6.8 ± 0.9 

^ Sum of the two Nal detectors. 

^ The 0.04 — 0.3 MeV count correspond to the Nal detector, while the others a sum of the 

Nal and Csl detectors. 

The flux is calculated by the data of Detector-B. 

Table 2. The obtained spectral parameters and x^/d.o.f. of the 071213. 

P. L. ^ Exp. P.L. ^ Exp. P.L. (fix) ^ 



a (MeV-^sr^^) 


(1.25 ±0.03) X 10^^ 


(7.19 ±0.02) X 10^" 


(4.9 ±0.6) X 10^ 




2.03 ±0.02 


1.88 ±0.01 


1.2 ±0.1 


ec (MeV) 
xLJd.o.f.'^ 




50 ±20 


7 


46.5/28 (1.6%) 


46.0/27 (1.3%) 


49.6/28 (0.72%) 


d (m) 


400li?° 


oc:q+210 
<JJU_240 


1501?° 



^ Power law model. 

^ Exponentially cut-off power law model. 

^ Exponentially cut-off power law with e^. being fixed at 7 MeV. 

Values in parentheses represent the probability with the given x^j^ and d.o.f.. 

^ Quoted errors of d are 90% confidence values, while other errors are 68% ones. 

Table 3. The obtained spectral parameters and x^/d.o.f. of 081225. 

P. L. Exp. P.L. Exp. P.L. (fix) 

a (MeV^^sr-i) (1.21 ± 0.09) x 10^ (6.4 ± 1.0) x 10^ (5.1 ± 0.4) x 10** 
(5 1.61 ±0.03 1.48 ±0.09 0.87 ±0.02 

ec (MeV) - 70 ± 80 7 

xl^Jd.o.f. 40.9/20 (0.38%) 40.6/19 (0.25%) 43.8/20 (0.16%) 

d{m) 300llg 250tii] 100+^^° ^ 

^ The lower value were unable to be determined. 

Table 4. The estimated for 071213 and 081225. 

071213 081225 

^max(deg.) if = 300 m H = 1000 m H = 300 m H = 1000 m 
15 1.9 X 10^° 6.2 X 10^ 3.7 x 10^ 1.2 x 10^ 
30 4.7 x 10^1 1.5 X 10" 9.3 x 10^° 3.0 x 10^° 
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